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Two-component Fermi and Bose gases with infinitely large interspecies s-wave scattering length 
o 3 exhibit a variety of intriguing properties. Among these are the scale invariance of two-component 
Fermi gases with equal masses, and the favorable scaling of Efimov features for two-component 
Bose gases and Bose-Fermi mixtures with unequal masses. This paper builds on our earlier work 
[D. Blume and K. M. Daily, arXiv: 1006.5002 and presents a detailed discussion of our studies 
of small unequal-mass two-component systems with infinite a 3 in the regime where three-body 
Efimov physics is absent. We report on non-universal few-body resonances. Just like with two- 
body systems on resonance, few-body systems have a zero-energy bound state in free space and a 
diverging generalized scattering length. Our calculations are performed within a non-perturbative 
microscopic framework and investigate the energetics and structural properties of small unequal- 
mass two-component systems as functions of the mass ratio k, and the numbers Ni and N2 of 
heavy and light atoms. For purely attractive Gaussian two-body interactions, we find that the 
(Ni,N2) = (2,1) and (3,1) systems exhibit three-body and four-body resonances at mass ratios 
k — 12.314(2) and 10.4(2), respectively. The three- and four-particle systems on resonance are 
found to be large. This suggests that the corresponding wave function has relatively small overlap 
with deeply-bound dimers, trimers or larger clusters and that the three- and four-body systems on 
resonance have a comparatively long lifetime. Thus, it seems feasible that the features discussed in 
this paper can be probed experimentally with present-day technology. 

PACS numbers: 03.75.Ss,05.30.Fk,34.50.-s 



I. INTRODUCTION 

Resonances arise in many different branches of physics. 
Near a resonance, certain physical observables are 
strongly enhanced, allowing one to probe the underlying 
physics at its "extreme". Resonances come in various 
types such as tunneling resonances (see, e.g., Ref. [l|), 
parametric resonances (see, e.g., Ref. [2}), stochastic res- 
onances (see, e.g., Ref. |3[), or Fano- Feshbach resonances 
(see, e.g., Ref. [4|]). Two-body Fano-Feshbach resonances 
play a key role in the study of ultracold atomic gases 
as they allow for the tuning of the interspecies and in- 
traspecies s-wave scattering length a s essentially at will 
while leaving all other system parameters essentially un- 
changed p-Q- In the ultracold regime, the scattering 
lengths determine the effective interaction strengths of 
the system. Thus, the tunability of a s allows for the re- 
alization of effectively repulsive and effectively attractive 
systems. Furthermore, the interaction strength can be 
tuned such that the system is effectively non-interacting 
or infinitely strongly interacting. Recent experiments 
that rely on this tunability include the study of the 
BCS-BEC crossover problem [9l4l5|. the polaron prob- 
lem [16[ , and the creation of bosonic and fermionic Fes- 
hbach molecules [l2|, These Feshbach molecules 
form an ideal starting point for creating, polar molecules 
in the rovibrational ground state |2l| - [23| . 

Fano-Feshbach resonances have also played a key 
role in recent experiments aimed at studying few-body 
physics. The most prominent example is the observation 
of Efimov physics [2J-|3l| , which relies on tuning the s- 



wave scattering length over several orders of magnitude. 
Efimov physics manifests itself most clearly in the vicin- 
ity of resonances. In particular, Efimov physics in the 
bosonic three-particle sector has been studied through 
the observation of enhanced three-body recombination 
loss rates near atom-dimer and atom-atom-atom reso- 
nances [U [H| • An atom-dimer resonance exists when 
the binding energy of the Efimov trimer equals that of 
the dimer. An atom-atom-atom resonance exists when 
the Efimov trimer has a zero-energy binding energy, 
i.e., sits right at threshold. Similarly, in the bosonic 
four-particle sector, the expected scaling associated with 
Efimov physics has been confirmed through measure- 
ments of enhanced four-body recombination loss rates 
at atom-atom-atom-atom, atom-trimer and dimer-dimer 
resonances [29l - l33| . An intriguing aspect of the features 
associated with Efimov physics is that the characteriza- 
tion of the three- and four-body sectors requires, for one- 
component Bose gases, just two parameters, the s-wave 
scattering length a s and a three-body parameter, which 
can for example be parametrized in terms of the binding 
energy of one of the Efimov trimers [H, |34| . These two 
parameters determine the positions of all universal reso- 
nances eluded to above. Efimov physics can also domi- 
nate the behavior of fermionic systems consisting of two 
or more components and of Bose-Fermi mixtures if the 
system parameters are tuned appropriately [35T - l4lj ] . 

This paper investigates few-body resonances for two- 
component systems with infinitely large interspecies s- 
wave scattering length a s in the regime where three-b ody 
Efimov physics is absent. Building on earlier work [3 a - 
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1371. |42h49| | . we address the following questions: (Ql) Un- 
der which conditions do atom-atom-atom resonances oc- 
cur for the (2, 1) system? (Q2) Under which conditions 
do atom-atom-atom-atom resonances occur for the (3, 1) 
system? (Q3) If the resonances discussed in Ql and Q2 
exist, what are their characteristics? (Q4) What are the 
differences and commonalities of system properties de- 
rived for finite-range (FR) and zero-range (ZR) interac- 
tions? (Q5) What are the implications of our theoretical 
studies for experiment [50|? 

The remainder of this paper is organized as fol- 
lows. Section III Al introduces the system Hamiltonian 
while Sees. Ill 51 and HTCl discuss two complementary ap- 
proaches for solving the time-independent Schrddinger 
equation for small two-component systems with unequal 
masses. Section UlTl presents our results for the energetics 
and selected structural properties for various parameter 
combinations. Lastly, Section IIVI summarizes our main 
results and concludes. A discussion of the main ideas 
and results of our work can be found in Ref. 48]. The 
present paper elaborates on the theoretical framework 
and provides a more detailed discussion of the results 
and their implications. In addition, the present paper 
presents structural properties, detailed comparisons be- 
tween observables derived within the numerical and an- 
alytical frameworks, and results for the (3, 2) and (4, 1) 
systems. 



II. THEORETICAL BACKGROUND 

This section introduces and discusses the system 
Hamiltonian that underlies our studies, and the 
techniques employed to solve the time-independent 
Schrodinger equation associated with this Hamiltonian. 
Our solutions to the Schrodinger equation are obtained 
following two distinctly different approaches. On the one 
hand, we pursue an analytical treatment that employs 
hyperspherical coordinates and determines the solution 
in terms of some unknown. On the other hand, we solve 
the Schrddinger equation numerically. The numerical so- 
lutions are then used to determine some of the unknowns 
that arise in the first approach, thereby providing an in- 
terpretation of the numerical results and, conversely, a 
check of the analytical framework. 



A. System Hamiltonian 

We consider two-component systems with N particles, 
where N = N% + under external spherically sym- 
metric harmonic confinement with angular trapping fre- 
quency uj. If Ni = 1 (i = 1 or 2), no permutation sym- 
metry of the ith species needs to be imposed; this im- 
plies that the single particle species can be fermionic or 
bosonic. In contrast, if Ni is greater than 1, the results 
depend on the permutation symmetry. Throughout, we 
impose fermionic symmetry for systems with N{ > 1. 



The masses of the two species are denoted by mi and 
777*2 ) and our model Hamiltonian H reads 

3=1 v 7 

£ G^ + ^?) + E E ^,).(1) 



j = l k=N x + l 



Here, fj denotes the position vector of the jth particle 
measured with respect to the trap center and Vth{fjk) 
with Tjk — \?j — rfe| the interspecies interaction poten- 
tial. Throughout, we assume that the particles of the 
same species are effectively non-interacting, i.e., we ne- 
glect intraspecies interactions. In the case of fermionic 
species this assumption is fullfilled with high accuracy 
for most systems since s-wave interactions are forbidden 
by symmetry and p-wave interactions are naturally sup- 
pressed [51 1 . 

Although nature provides us with only a finite number 
of discrete mass ratios k, where 



K = ?77l/7772, 



(2) 



we find it useful to treat k as a continuous variable 
to unravel the key physics, i.e., to see how the physics 
changes as a particularly interesting mass ratio is ap- 
proached. Experimentally, the effective mass ratio of 
two-component systems could be tuned by loading the 
system into an optical lattice [13]. In addition to the 
mass ratio, we vary the number of heavy and light parti- 
cles, the angular momentum and parity of the state under 
consideration, and the two-body interaction Vth- 

Our calculations are performed for two classes of in- 
teraction potentials, FR and ZR interactions. Our FR 
calculations are performed for a purely attractive Gaus- 
sian potential V g with range rg and depth Vq (Vq > 0), 



v s( r ) = -Voexp 



V2r 



(3) 



While the majority of our FR calculations considers the 
infinite scattering length limit, i.e., l/a s = 0, we also 
explore how the system behavior changes as \a s \ is de- 
creased. In particular, we pick ro and then adjust the 
depth Vq so that the two-body potential has the desired 
free-space s-wave scattering length a s . For positive a s , 
we restrict ourselves to potentials that support a single 
s-wave two-body bound state. For negative a s , we re- 
strict ourselves to potentials that support no two-body 
bound state. Our ZR calculations employ the Fermi- 
Huang pseudopotential U zr [52|, HU , 



2nh 2 a s d 

U zr (r) = 5{r)—r, 

fi or 

where /i denotes the reduced mass, 

/! = 777 1 7772/ (mi + 7772). 



(4) 



(5) 
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Our ZR calculations are restricted to unitarity, i.e., to 
the regime where the s-wave scattering length is infinitely 
large and where the two-body system in free space is at 
the verge of supporting a zero-energy bound state. 

The system under study is characterized by the follow- 
ing length scales: the harmonic oscillator length aho, 



dho = \A/(2/xw), 



(6) 



the range rg of the interaction potential (r = for V ZI ), 
and the s-wave scattering length a s . At unitarity, the 
s-wave scattering length no longer defines a meaningful 
length scale. For sufficiently large k, an additional length 
scale is given by the generalized ./V-body scattering length 
(see Sec. Ill B| ). Our FR calculations are performed in the 
regime where ro <C ai ao . In fact, one of the goals of this 
paper is to quantify how observables obtained for FR 
interactions approach those determined in the ZR limit. 
As will be shown in Sec. lIIIl FR effects can be appreciable 
for two-component unequal-mass systems. 

Our approaches outlined in Sees. Ill Bl and III CI take 
advantage of the fact that the center-of-mass degrees of 
freedom i? cm separate off. Throughout, we assume that 
the center-of-mass wave function is in the ground state, 
and we label the solutions *5> to the time-independent 
Schrodinger equation in the relative coordinates by the 
relative orbital angular momentum L and the relative 
parity II. The corresponding relative eigenenergies will 
be denoted by E. The Hamiltonian given in Eq. ([T]) de- 
scribes the trapped system. As we discuss in Sec. Ill Bl the 
eigenenergies of the trapped iV-body system can be ex- 
pressed in terms of the generalized scattering length that 
characterizes the corresponding iV-body scattering prob- 
lem in free space. This connection between the trapped 
and free-space systems, which is of course well known 
for the two-body problem [54j j . provides a great deal of 
insight. 



B. Hyperspherical coordinate treatment 



Within the hyperspherical framework [55 58], the 
3N — 3 relative coordinates are divided into 3N — 4 hy- 
perangles, collectively denoted by J7, and a single length, 

the hyperradius R, /iR 2 = Y^jLi m j(^j ~ Rem) 2 - In 
the present context, these hyperspherical coordinates 
are particularly appealing since the relative wave func- 
tion of the unitary system with ZR interactions has 
been shown to separate for any number of particles 
into a hyperradial part F uq (R) and a hyperangular part 
[42], V uq (R,n) = R-^ N -^' 2 F vq {R)^ v {n). Here, 
v and q denote hyperangular and hyperradial quantum 
numbers for a given angular momentum L and parity 
II; v takes the values 0, 1, ■ • • while q takes non-integer 
values (see below). For ZR interactions with l/a s = 0, 
the relative Schrodinger equation can therefore be solved 
in a two-step process. First, the hyperangular functions 
$1,(0), or so-called channel functions, are determined 




FIG. 1: (Color online) Effective hyperradial potential curve 
V Vt eff(R), Eq. ©, as a function of R for, from top to bottom, 
s„ = 1,3/4, 1/2,9/20 and 0. 



by solving the hyperangular Schrodinger equation [42j. 
The corresponding eigenvalues are related to the coeffi- 
cients s„, which determine the effective potential curves 



2fiR 2 



-^ 2 R 2 . 



(7) 



Figure Q] shows the effective potential curves V v , e ff(R) for 
different s v . For s„ = 1/2, the hyperangular 1/R 2 term 
vanishes and V„ je ff (i?) reduces to the trapping potential, 
i.e., the second term on the right hand side of Eq. (J7J. 
For s v > 1/2 and s„ < 1/2, the hyperangular 1/R 2 term 
is repulsive and attractive, respectively, and dominates at 
small R. Second, the hyperradial Schrodinger equation 



+ VuMR) F »«( R ) = E„ q F uq {R) (8) 



-h 2 d 2 

2fi OR 2 

is solved for F vq (R) and E vq . If one uses FR instead of 
ZR interactions, the hyperangular and hyperradial parts 
of the wave function vp, in general, do not fully sepa- 
rate, implying non-vanishing coupling matrix elements 
between the different hyperangular channel functions. In 
a first approximation, however, these couplings can be 
neglected if ro -C &ho- The framework outlined in this 
section is thus not only applicable to ZR interaction but 
also, at least within an approximative scheme, to FR in- 
teractions (see Refs. [59l. t60T] ) . 

Solving the hyperangular Schrodinger equation, i.e., 
determining the s„, is, in general, a non-trivial task. 
For the three-body system with ZR interactions and in- 
finitely large a s , however, the eigenvalues can be obtained 
by solving a simple transcendental equation for each L n 
symmetry 0, [2I, |H, S3, [Hj]. Figure M shows the s 
coefficients for the (2,1) system, i.e., the system with 
two heavy fermions and one light atom, at unitarity with 
L = — 3 and II = (— 1) L as a function of the mass ra- 
tio k. The state with the smallest sq value has L u = 1~ 
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FIG. 2: (Color online) so coefficient as a function of k for the 
(2, 1) system at unitarity with ZR interactions and L n = + 
(dotted line), L n = I - (solid line), L n = 2+ (dashed line) 
and L n = 3~ (dash-dotted line). Three-body Efimov physics 
is absent if so > 0, which corresponds to k < 13.607 for the 
L n = 1~ system and to k < 76.0 for the L n = 3~ system. 



symmetry for mass ratios k < 13.607 (solid line in Fig. [5]). 
The coefficient s$ decreases from 1.773 for k = 1 to for 
k ss 13.607. For later reference, we note that Sq equals 1 
for k « 8.619 and 1/2 for k = 12.313. For k > 13.607, 
so becomes purel y imag inary and Efimov physics comes 
into play [111, HH, l35l437j . The sq coefficient of other odd 
L states (the dash-dotted line shows the so coefficient 
for the L n = 3~ state) becomes imaginary for much 
larger mass ratios. For even L states, s increases with 
increasing k and approaches an integer value in the large 
k limit. For larger systems with infinitely large inter- 
species s-wave scattering length, the s v coefficients are 
not known in general. 

We now consider the hyperradial Schrodinger equation, 
Eq. J5]), for the TV-body system, which can be solved for 
arbitrary s v . Generalizing the quantum defect theory 
type approach from Ref. [62| from s„ = 1/2 to arbitrary 
s v , the solution F uq (R) to the second order differential 
equation can be written as 



F vq (R) = N vq [f vg (R) - ta,n(TTfi ug )g„q(R)} 



(9) 



where N vq denotes a normalization constant. The 
quantum defect [L vq determines the relative contribu- 
tions of the regular solution f vq and the irregular so- 
lution g vq , f vq {x) = A vq x s » +1/2 exp(-x 2 /2) 1 F 1 (-q,s v + 
l,x 2 ) and g uq {x) = B vq x~ s " +1 / 2 exp(— x 2 /2)iFi(— q_— 
Su, — s u + l,x 2 ). Here, A uq and B uq denote constants [63J, 
and x the dimensionless hyperradial coordinate, x = 
R/(V2a ho ). The non-integer quantum number q is re- 
lated to the eigenenergy E vq through 



E vq = {2q + s v + l)fvjJ. 



(10) 



To determine the allowed values of (j, uq , we enforce that 
F vq {x) vanishes at large x, resulting in the condition 



sin[7r(/i„ g + q)] — 0. This condition allows fi uq to be 
eliminated, leaving the quantum number q as the only 
unknown. 

The quantization condition, i.e., the allowed q values, 
are determined by investigating the small x behavior of 
F vq {x). For x — > 0, f vq {x) behaves as x s,y+1 ^ 2 and is 
well-behaved or less strongly diverging than x~ x / 2 for 
s v > —1 while g vq (x) behaves as x ~ Sv+1 / 2 and diverges 
faster than x~ x / 2 for s y > 1. Thus, we consider the 
regimes s u > 1 and < s„ < 1 separately (the Efimov 
regime with imaginary s v is not treated in this paper). 
For s v > 1, g V q{x) must be eliminated [64| and F uq (x) is 
determined by the exponentially decaying piece of f uq {x). 
The quantization condition becomes, in agreement with 
Ref. [42], q = 0, 1, • • • and the corresponding energy is 



referred to as E 



E i>vq = (2q + s v + i)hw; q = 0,l, 



(11) 



For < s v < 1, both f vq {x) and g vq {x) are well-behaved 
and the allowed q values depend on the boundary condi- 
tion at small x. This boundary condition is determined 
by the true atom-atom interactions and cannot be de- 
rived within the ZR framework. Similarly to the case 
of Efimov trimers [24T - I261 ] . the value of the short-range 
hyperradial boundary condition can be thought of as an 
extra parameter that is needed to specify the solution. In 
our case, this parameter characterizes the iV-body system 
and can be parameterized, e.g., by the logarithmic deriva- 
tive L vq {x Q ), L uq (xu) = [{dF uq {x) / dx) / F uq (x)] x=Xo . 

For s„ > (s„ not equal to an integer), the normal- 
ized hyperradial wave function F vq (R) can be compactly 
written as 



N vq exp 



F uq {x) = 
^+i/2 U (-q,s v + l,x 2 ), (12) 



where U denotes the confluent hypergeometric function 
of the second kind and 



I -2 S m(Trs v )T{l-q)T(-q-s v ) 
7r + ir 2 qcot(irq) + qir[i/)(q) — tp{—q ~ s u )] 



;(i3) 



here, ip denotes the digamma function [651 ] . Using 
Eq. p2[) . the logarithmic derivative at x — xq can be 
written as 



L vq (x ) 



x 



2{q + s v ) 



x 

U(-q,s v ,x%) 



xq 



U(-q,s u + l,xi) 



(14) 



Figure Ufa) shows the quantum number q for s v = 3/5 
for x < 1 (i.e., x a = 1 x 10 -3 , 2 x 10~ 3 and 3 x 10 -3 ) as 
a function of the inverse of the logarithmic derivative, q 
is approximately equal to integer values (q s=s 0, 1, • • ■ ) for 
nearly all values of the logarithmic derivative. Exceptions 
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the ZR limit, in agreement with note [43] of Ref. [42 
We refer to the corresponding energies as E g , vq , 



cr 



-1 - 



\ 
V 


1 


•. s 

\ 
V 


- 


-. -s" 1 - 

1 

: 1 

■ (a) \ 1 
. ■ . i - . ri 





Eg, vq — (2q 



l)fiw; q = —s„, —s v + 1, • 



(15) 



105 -0.04 -0.03 -0.02 -0.01 



cr 0- 



-1 - 



i i i i K^J 














-(b) 
■ i i 


, 1 i , , 





1 

<y 
-1 



'U-f-..'. - i — i — r 

\ *'••. 



■(c) 



\ 
\ 
I 
I 

J_L 



0.01 0.02 0.03 0.04 0.05 
1/L (x n ) 

FIG. 3: (Color online) Non-integer quantum number q as a 
function of [L vq (xo)]~ for (a) s v = 3/5, (b) s v — 1/2, and 
(c) s v — 2/5. In panels (a) and (c), the solid, dashed and 
dotted lines correspond to io = 1 x 10 -3 , 2 x 10 -3 and 3 x 
10 -3 , respectively. In panel (b), the solid, dashed and dotted 
lines correspond to xo = 10 -5 , 10 -3 and 10 _1 , respectively; 
the solid and dashed lines are indistinguishable on the scale 
shown. Note that the scale of the :r-axis is different in all 
three panels. 



occur for negative [L vq {xQ)]^ 1 values with \\/L vq {xo)\ 
small, where q drops by roughly 1 over a small range of 
[L vq (xo)]^ 1 values. The value of \l/L„ q (xo)\ at which 
q drops decreases as Xq decreases and reaches in the 
ZR limit. The deviation of q from w 0, 1, ■ ■ • sig- 
nals the existence of an TV-body resonance, which arises 
when F vq (x) is fully determined by the irregular solution 
g vq {x). When only g vq (x) contributes, we find that the 
quantization condition becomes q = — s v , — s u + 1, • • ■ in 



For s„ = 3/5, q equals —3/5, 2/5, • • • on resonance and 
\1/L vq (x )\ = for x Q -> 0. 

The dependence of q on the inverse logarithmic deriva- 
tive l/L vq (xo) for s v = 2/5 [Fig. EJc)] is similar to 
that for s v = 3/5 [Fig. E^a)], with the main difference 
that q now drops sharply for small positive l/L vq {xo) 
as opposed to negative l/L vq (xo). For s v = 2/5, q 
equals — 2/5, 3/5,- •• on resonance and \l/L vq (xo)\ = 
for xq — > 0. We note that the trapped iV-body sys- 
tem supports a deep-lying bound state if the logarithmic 
derivative is negative [not shown in Fig. [3fc)]. 

The dependence of q on [L^ q (xo)}^ 1 for s v = 1/2 
[Fig. Hfb)] is distinctly different from that for s„ = 3/5 
and 2/5. Figure [3] shows that the q values depend much 
more weakly on xo for a v = 1/2 than for s v ^ 1/2. Fur- 
thermore, q approximately equals 0, 1, • • • for 1/L u (x ) = 
0. As \1/L uq (x )\ increases, the q values change gradu- 
ally and approach half-integer values for large | X/L vq (xq ) \ 
(on resonance, q = — 1/2, 1/2, • • • ). In summary, for 
s„ 7^ 1/2, the system exhibits an iV-body resonance in 
the ZR limit for l/L uq (xo) = 0; for q = —s v , e.g., the 
third term on the right hand side of Eq. (f!4|) vanishes 
and lim Xo _>.o |L„ g (xo)| = oo. For s v = 1/2, in contrast, 
the system exhibits an ./V-body resonance in the ZR limit 
for l/L uq (xo) = oo; for q = —s u , e.g., the first and the 
third term on the right hand side of Eq. (IT41 vanish and 
lim Xo _>.o L ug (x ) = 0. 

We now introduce a framework that expresses the 
resonance condition in terms of the generalized energy- 
dependent scattering length V Sv (E). To define V Sl/ (E), 
we consider the free-space system, i.e., we set the trap- 
ping frequency uj to 0. Compared to the trapped sys- 
tem, only the hyperradial Schrodinger equation changes 
[see Eqs. J?]) and ©]. The hyperradial solution for the 
free-space system with positive energy E can be written 
analogously to the solution for the trapped system [see 
Eq. ©], 

F vk {R) = 

N vk [r 1/2 J Sv (kR) - tan(5 s „ {k))R 1/2 Y Sv (feii)] , (16) 

where we use the continuous variable k to label the solu- 
tion of the hyperradial Schrodinger equation; k is defined 
in terms of the iV-body s cattering energy E and the hy- 
perradial mass n, k = y / 2fiE/h 2 . In Eq. ([To]) . J Sv and 
Y Sv denote the Bessel functions of the first and second 
kind, respectively. The phase shift 5 Su (fc) characterizes 
the iV-body scattering process and can be used to de- 
fine the generalized energy-dependent TV-body scattering 
length V Sv (k), 



VsAk) 



tan(^(fc))2 2s -r(s t/ )r(a t/ + l) i 

k 2s " 7T 



(17) 
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throughout this paper, the generalized energy-dependent 
scattering length is written, depending on the context, 
as V Sv (k) or V S „(E). The power of k in the denominator 
on the right hand side of Eq. (TTTj) is chosen such that 
the generalized energy-independent iV-body scattering 
length V Sv (0), V s „(0) = linifc^o V Sv (k), is well-behaved, 
i.e., such that V s „(k) approaches a constant in the k — > 
limit. The "extra" factors on the right hand side of 
Eq. (1TH) are chosen such that V Sl/ (0) reduces to R 2 ?" if 
the boundary condition F„k{Ro) = is imposed. We note 
that the definition of V Sv (E), Eq. ([X7]l . has some similari- 
ties with that employed to describe the scattering of two 
particles with vanishing azimuthal quantum number in 
two spatial dimensions [Hil For s v — 1/2, the N- 
body hyperradial scattering problem becomes formally 
identical to the scattering between two three-dimensional 
s-wave interacting particles and V Sv (0) becomes formally 
equivalent to the usual three-dimensional s-wave scat- 
tering length a s . In general, however, the generalized 
scattering length has units of (length) 2s " . 

To express the resonance condition for the trapped sys- 
tem in terms of the generalized energy- dependent scat- 
tering length V S „(E), we relate the eigenenergies of the 
trapped A^-body system to V s „ (E). To this end, we calcu- 
late the logarithmic derivative of the free-space solution 
F v k{R), Eq. (fl"6]) . for small xq. Considering < s v < 1 
and keeping terms up to order x 2s "~ 1 , we find 



xo V Su {E) 
2- s - +1 7rcot(7rs l ,) I E 
foui 



(18) 



For s„ = 1/2, the first and third terms on the right 
hand side of Eq. (|18p vanish and we find L u k(xo) rs 
~V2a ho /V Sv (E) or L vk (Ro) « -[V s „ (E)}-\ To ob- 
tain an explicit relationship between the eigenenergies of 
the trapped iV-body system and the generalized energy- 
dependent scattering length V Sv {E) for < s v < 1, we 
expand the logarithmic derivative of the trapped system, 
Eq. flT4|) . up to order Xq S " _1 and set it equal to the log- 
arithmic derivative of the free-space system [right hand 
side of Eq. (Tig])]. 



2s v 
l ho 



VsAE) 
where 



+ h(E,s u ) 



r( 



h(E,s u ) 



7rcot(7rs t/ ) / E 

2 2s -r( Sl/ )r(i + S „) \nu 



(19) 



(20) 



It can be seen that the 1/xo divergencies are can- 
celed and, furthermore, that Eq. (TH)|) is — at this level 
of approximation — independent of xq. Equation (fig)) 
determines the energy of the trapped iV-body system 
in terms of the generalized energy- dependent scattering 
length V Sv (E) and can be solved self-consistently for the 
eigenenergies E. 



S3" 




FIG. 4: (Color online) Solid, dashed, dotted and dash-dotted 
lines show the quantity h(E, s v ), Eq. (|20[) . as a function of s„ 
for E = 1 x 10~ 2 huj, 1 x 10 _1 fiw, lhw, and lOftu, respectively. 



For s v = 1/2, h(E,s„) vanishes and T(l - s„)/[T(l + 
s„)2 s "] reduces to y/2. Setting [L„ q (xo)] _1 equal 
to — V Sv (i?)/(v / 2aho) [see discussion after Eq. ([15])]. 
Fig. GUb) can be interpreted as showing the non-integer 
quantum number q of the trapped -/V-body system as a 
function of — V s „(-E)/(v2ciho)- A vanishing generalized 
energy-dependent scattering length implies an infinitely 
large logarithmic derivative, and an infinitely large gener- 
alized energy-dependent scattering length implies a van- 
ishing logarithmic derivative. From Eq. (|19p . it follows 
for s„ = 1/2 that E = E g>l/q if V S „(E) diverges, i.e., 
the Af-body resonance occurs when [V Siy (-E)] _1 vanishes. 
The analysis outlined here for s„ = 1/2 is formally iden- 
tical to that of the trapped two-particle system. Identify- 
ing V Sv (E) with the s-wave atom-atom scattering length, 
Eq. (|19p is identical to the well-known eigenequation for 
two trapped s-wave interacting particles with ZR inter- 
actions |54j . 

For s v ^ 1/2 and E > 0, h(E, s„) does not vanish and 
introduces an additional energy dependence on the left 
hand side of Eq. (|19p . which originates from the explicit 
energy-dependence of the logarithmic derivative of the 
free-space solution. If we artificially set h(E, s„) to zero, 
the eigenenergy E equals E g>vq when V Su (E) diverges. 
Inclusion of h(E, s„) shifts the energy for all V Sv (E) ^ 
(s„ =/= 1/2) and modifies the resonance condition. In 
particular, the resonance condition becomes 



(21) 



Figure @] shows h(E,s v ) as a function of s u for four dif- 
ferent energies, i.e., E = 10~ 2 hu — 10hu. h(E,s v ) van- 
ishes if Sy equals 1/2 (see also above), and increases with 
increasing \s v — 1/2|. For E > 0, h(E, s v ) takes on neg- 
ative values for 1/2 < s„ < 1 and positive values for 
< s u < 1/2. As can be seen, h(E,s v ) approaches 
— oo as s v approaches 1. Thus, for large |V a „ {E)/a 2 ^ \ 
and s v 7^ 1/2, the ar*" /V Su (E) term is small and the 
h(E,Sy) term dominates the left hand side of Eq. (|19l) . 
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FIG. 5: (Color online) Dashed lines show the non-integer 
quantum number q of the trapped iV-body system, obtained 
by solving Eq. (|19p self-consistently, as a function of the gen- 
eralized energy- dependent scattering length V s „ {E) for (a) 
s v = 3/5 and (b) s„ — 2/5. For comparison, solid lines show 
the g values that result when h(E, s v ) is artificially set to zero. 



ters of a lattice confinement 44]. The realization of an 
TV-body resonance for s v ^ 1/2 requires extreme fine- 
tuning since both h(E, s„) and V S „(E) vary with energy. 
On the other hand, for s„ = 1/2, only V s „ varies with 
energy. The energy-dependence of V 8v (E) is expected to 
be weakest for s u = 1/2 since the effective angular mo- 
mentum barrier in the R coordinate vanishes in this case. 
This reasoning is motivated by what is known from "ordi- 
nary" two-particle scattering (see, e.g., Ref. [HI): In gen- 
eral, two-particle s-wave scattering (no angular momen- 
tum barrier) exhibits a much weaker energy-dependence 
than two-particle p- or d-wave scattering (finite angular 
momentum barrier). We conclude that it should be more 
likely to realize an iV-body resonance if s„ = 1/2 than if 
s v ? 1/2. 

Lastly, we note that application of the two-body scat- 
tering framework to the hyperradial problem at hand 
implies that the free-space system supports an iV-body 
bound state with zero energy when V Sv (0) diverges. Fur- 
thermore, the iV-body system supports a single weakly- 
bound free-space bound state when V Su (E) is large and 
positive, and no weakly-bound free-space bound state 
when V Su (E) is negative. 

The following section introduces the stochastic varia- 
tional (SV) approach which we use to numerically solve 
the Schrodinger equation in the relative coordinates for 
few-body systems interacting through FR potentials. 
This approach employs Jacobi coordinates instead of hy- 
perspherical coordinates and makes no assumption about 
the small R behavior of F uq (R): the small R behavior is 
not treated as input but instead is a natural part of the 
solution. Connections between the results based on the 
approaches discussed in Sees. Ill Bl and III Cl will be made 
in Sec. HIH 



For small | V a „ (E)/a^ | , in contrast, the a^ " /V S „(E) 
term is large and the h(E, s v ) term on the left hand side 
of Eq. (fl"9| can be neglected. Figures Ufa) and (b) show 
the eigenenergy of the trapped system as a function of the 
generalized energy-dependent scattering length V Su {E) 
for s„ = 3/5 and 2/5, respectively. Dashed lines show 
the non-integer quantum number q obtained by solving 
Eq. (IT51) self-consistently while solid lines show the non- 
integer quantum number q that results when h(E, s u ) is 
artificially set to 0. Figure [5] shows that the inclusion of 
the h(E, s v ) term shifts the eigenenergies up (down) com- 
pared to those calculated for h(E,s„) = for s v = 3/5 
(s„ = 2/5). The term h(E,s v ) has a small effect when 
|V S „(-E)| is small but introduces a notable shift of the 
energies when |V S „(-E)| is large. 

As discussed after Eq. (jTTJ) , the ZR framework em- 
ployed in this section does not allow for the determi- 
nation of the actual value of the logarithmic derivative 
or, equivalently, the generalized scattering length. Let 
us imagine that the generalized scattering length V s „ (E) 
can be controlled experimentally by tuning the interac- 
tions between like atoms [44| or by varying the parame- 



C. Stochastic variational treatment 

The SV approach li^ Ff(I expands the relative wave 
function ^ in terms of basis functions ifk , 

*= c k A[ip k (x)}, (22) 
fe=i 

where the Ck denote expansion coefficients, A denotes 
an operator that ensures the anti-symmetry of the basis 
functions, and x collectively denotes the Jacobi vectors 
pj, where j = l,--- ,N — 1. The basis functions [see 
Eqs. (|2"3")l and (|2"5j) for their explicit forms] are chosen such 
that the Hamiltonian matrix can be constructed analyt- 
ically. The eigenenergies E of the relative Hamiltonian, 
which provide an upper bound to the exact eigenener- 
gies, are then obtained by diagonalizing the Hamiltonian 
matrix. The variational bound can be improved system- 
atically by increasing the size of the basis set (i.e., by in- 
creasing -/Vbasis) and by varying the parameters and 
d\f [see Eq. (g31)] or s< fc ) and d\f [see Eq. p5])]. Our pri- 
mary interest in this work is to describe the energetically 
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lowest-lying state of various (Aq, N2) systems. The ener- 
getically lowest lying state of the (2, 1) system, e.g., has 
natural parity (see Fig. [2J while that of the (3, 1) system 
has unnatural parity. Thus, we need to consider basis 
functions that can describe both natural and unnatural 
parity states. 

To treat natural parity states, i.e., states with II = 
(— we employ basis functions that are written in 
terms of the spherical harmonic Ylo, which determines 
the relative orbital angular momentum of the system and 
depends on the unit vector 68], 



tp k (x) = \v(V\ L Y L0 (vW) 



exp 



N 

■E 



i<3 \^ d 'ij 



(k) 



where 



JV-l 

i=i 



,(23) 



(24) 



Here, the Uj form a (JV— l)-dimensional parameter vec- 
tor that determines how the relative orbital angular mo- 
mentum L of the {N\, N2) system is distributed among 
the (N — 1) Jacobi vectors pj. We find that the optimal 

set of widths d\) depends quite strongly on the mass ra- 
tio k. For the (2, 1) system with equal-masses, e.g., three- 
body bound states are absent [H, [7l| • This implies that 
three-body correlations are largely absent, and that the 

(k) 

contribution of basis functions with more than one d)- of 

l j 

the order of the range of the underlying two-body poten- 

tial contribute negligibly to the wave function [gO, [z3 ■ 

When k as 12, in contrast, three-body correlations are 

non-negligible and basis functions ipk that are character- 
ful 

ized by three widths d\j of the order of the range of the 
two-body potential contribute notably |48| . 

To describe states with unnatural parity, we employ so- 
called geminal-type basis functions tfk that are neither 
eigenfunctions of the angular momentum operator nor 
the parity operator [68|, , 



(fik(x) = exp 




(25) 



Here, (s«) T x is just the dot product between two 
3(N — 1) dimensional vectors. The 3(N — l)-dimensional 
parameter vector s^^ is, together with the N(N — l)/2 

widths djj , optimized semi-stochastically for each basis 
function ipk. In general, the determination of the natural 
parity states [see Eq. (|23| ] is, for the same N, numeri- 
cally signiheantly more efficient than that of unnatural 
parity states [see Eq. (|25])]. 

The SV approach also allows for the determination of 
structural properties. To calculate structural properties, 
we follow two different approaches [73[ : (i) We implement 
the analytically known matrix elements [68j for the op- 
erator A of interest, and calculate the expectation value 



of A using the known c^, where k = , Abasis- (ii) 

We calculate the expectation value of the operator A by 
sampling the normalized density \ifj\ 2 through Metropo- 
lis sampling [74[. In the limit of infinitely many Monte 
Carlo samples, the results of approach (ii) should agree 
with those of approach (i). Where possible, we have used 
this to check our implementations. To calculate the hy- 
perradial density Phypcr(-R), where 



Phyper(R)dR — 1, 



(26) 



we employ approach (ii) . The pair distribution functions 
Ph\(f) and Phhir) for heavy- light and heavy-heavy atom 
pairs, normalized such that 



4-7T / Pw.(r)r 2 dr = 1 



(27) 



[and similarly for Phh (?")], is calculated using approach 
(i) for the wave function written in terms of the basis 
functions given in Eq. (|23[) and approach (ii) for that 
written in terms of the basis functions given in Eq. (|25l) . 



III. RESULTS 

This section presents our results from the SV calcu- 
lations for FR interactions and interprets our findings 
within the hyperspherical framework. Sections 1111 Al and 
IIIIB I present the energetics and structural properties for 
the (2, 1) and (3, 1) systems, respectively. 



A. Three-body resonances 

We start our discussion with the infinitely strongly in- 
teracting (2,1) system at unitarity, and then discuss the 
behavior away from the two-body resonance. Since the 
L n = 1~ state is, as discussed in Sec. Ill BL the energet- 
ically lowest lying state for all n of interest, we restrict 
our SV calculations to this symmetry. For each mass ra- 
tio k, we consider different tq with tq <C &ho- Symbols 
in Figs. [6] and [7] show examples for the ground state en- 
ergy of systems with k — 1 through 12.316 interacting 
through V g . By fitting the range-dependent SV energies 
to simple three- to five-parameter expressions, we extrap- 
olate the FR SV energies to the ro — > limit. For n = 1 
and 6 (circles and squares in Fig. \§§ , the SV energies ap- 
proach the ZR limit approximately linearly from above 
and below, respectively. For n = 10, the FR energies are 
best described by a quadratic three-parameter fit. As 
k increases further, the range-dependence increases no- 
tably (see circles and squares in Fig. [7] for n = 12.25 
and 12.3). For k = 12.3131 and 12.314 (diamonds and 
triangles in Fig. [7]), in contrast, the range-dependence is 
comparatively small and the SV energies approach the 
ZR limit approximately linearly from above. Finally, for 
k = 12.316 (crosses in Fig. [7]) the energies decrease as 
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FIG. 6: (Color online) Circles, squares and diamonds show the 
SV energies for the (2, 1) system with L n = 1~ at unitarity 
for k = 1, 6 and 10, respectively. The unlike particles interact 
through V g . Solid lines show three-parameter fits of the form 
Y^ij=o c i r o to the SV energies. 
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FIG. 8: (Color online) Energies of the (2, 1) system with 
L n — V at unitarity. Circles, crosses and pluses show the 
SV energies for V g , extrapolated to the ZR limit, for the three 
energetically lowest-lying states as a function of k. Solid and 
dotted lines show Ef,o q (q = and 1) and -E s ,o<j (q ~ —so 
and —so + 1), respectively. This figure has been adapted from 
Ref. HH]. 
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FIG. 7: (Color online) Circles, squares, diamonds, tri- 
angles, pluses and crosses show the SV energies for the 
(2, 1) system with L n — 1~ at unitarity for k — 
12.25,12.3,12.3131,12.314,12.315 and 12.316, respectively. 
The unlike particles interact through Vg. Solid lines show 
three- to five-parameter fits to the SV energies. 



l/r§ for sufficiently small ro, indicating the presence of 
a three-body bound state. For yet larger k (not shown 
in Fig. [7]), the energy of the trapped system becomes 
negative. The appearance of the three-body bound state 
indicates the presence of a three-body resonance. We es- 
timate the resonance for the potential V g in the ro — > 
limit to be located at k w 12.314(2). The uncertainty 
in our estimate for the resonance position arises mainly 
from the extrapolation of our SV energies to the ZR limit. 
The SV energies for n = 12.315, e.g., are nearly equally 
well described by fits of the form c_2/Vq + c_i/V + c 
(shown in Fig. [7J and cq + c\r$ + C2Tq (not shown in 
Fig.©. 



Symbols in Fig. [8] show the extrapolated ZR ener- 
gies for the three energetically lowest lying states of the 
(2, 1) system with L n = 1~ interacting through V g with 
l/a s — 0. The dropping of the energies around k « 12.3 
associated with the three-body resonance is clearly visi- 
ble: The ground state energy (circles) becomes negative 
while the energies of the first excited state (crosses) and 
the second excited state (pluses) drops by approximately 
2fku. For comparison, solid and dotted lines in Fig. [8] 
show the ZR energies Ef$ q [Eq. (fTT| with q = and 1] 
and Eg t Q q [Eq. (|15l) with q = —sq and — sq + 1], respec- 
tively. Away from the three-body resonance, the extrap- 
olated ZR energies agree well with Ef, 0q . On resonance, 
the extrapolated ZR energies agree well with E g .o q . Our 
analysis is fully consistent with the general discussion of 
Ref. [42j , where it was found that universal states for sys- 
tems under spherically symmetric harmonic confinement 
must have an energy larger than lhaj and that states 
with energy less than Ihui are necessarily non-universal. 
The admixture of the irregular solution g vq requires that 
the boundary condition of the hyperradial wave function 
be specified, which makes the system properties depen- 
dent on an additional parameter and thus non-universal. 
It has been pointed out previously [42j that the three- 
body system is again scale-invariant at the three-body 
resonance. Our work (see also Ref. (48|) shows an ex- 
ample for this non-universal three-body resonance. This 
non-trivial three-body resonance has also recently been 
investigated by Gandolfi and Carlson [4{|> who studied 
the free-space problem. 

The mass ratio at which the three-body resonance ex- 
ists depends on the details of the underlying two-body 
potential. As has been argued in Ref. [48| and in Sec. Ill Bl 
it seems most likely that the three-body resonance occurs 
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FIG. 9: (Color online) Symbols show the SV energies for the 
(2, 1) system with L n — 1~ interacting through V g at unitar- 
ity for (a) ft = 12 and (b) k — 12.7, respectively, as a function 
of rn. For k = 12, the energies of the ground state (circles) 
and first excited state (pluses) are shown. For k = 12.7, 
the energies of the first and second excited states (pluses and 
squares) are shown. Solid lines show fits to the SV energies. 
Horizontal lines show the energy -E/,oq for q = and 1. 



when the effective hyperangular momentum barrier van- 
ishes, i.e., when s„ w 1/2. For the (2,1) system with 
ZR interactions and L u — 1~ symmetry, this is the case 
when k w 12.3131 (see Fig. [5]); this value is close to the 
resonance position found for the purely attractive Gaus- 
sian potential V g . In our example, the trimer on reso- 
nance is large [this follows from Eq. ((12)) ]: this implies a 
long lifetime and thus opens the intriguing possibility of 
studying trimer correlations in many-body systems. 

Symbols in Figs. (Hla) and (b) shows the SV energies 
for the (2, 1) system with k — 12 and 12.7, respectively, 
interacting through V g with l/a s = as a function of 
ro. For both k, the SV energies for the "lowest universal 
state" and the "second lowest universal state" are shown 
[for k = 12.7 ', the ground state (not shown) corresponds 
to a tightly bound non- universal trimer]. The extrapo- 
lated ZR energies of the two universal states agree well 
with the energies £7,0? (<Z — arL d 1), which are shown 
by horizontal lines. This implies that the 2hu) spacing 
expected for universal states 42[ is fullfilled with good 



accuracy. For k = 12 and 12.7, our fits result in an en- 
ergy spacing of 2.014fiw and 1.9847kj, respectively. For 
finite ro, the spacing between the two lowest universal 
states deviates from the 2huj spacing. For ro ~ O.Olaho, 
e.g., we find a spacing of ps l.fJOfiu; and 2.13fiw for k = 12 
and 12.7, respectively. 

To make an explicit connection between the FR SV en- 
ergies and the energies obtained within the hyperspheri- 
cal framework (see Sec. Ill B| ). we assume a direct propor- 
tionality between the range ro of the two-body interac- 
tion potential and the hyperradius Ro at which the log- 
arithmic derivative of the hyperradial function F uq (R) is 
imposed. The proportionality factor between ro and Rq 
cannot be determined within the model; we find that a 
proportionality factor of order 5 to 10, i.e., r ~ 5 — 10i?o, 
is appropriate. Figures [R)T a)-(d) exemplarily show the 
eigenenergies predicted by the hyperspherical framework 
as a function of Ro for s v s» 0.5579 and the lowest allowed 
q value, for s v = 1/2 and the lowest allowed q value, and 
for s v — 0.4180 and the lowest and second lowest allowed 
q values for various values of the logarithmic derivative. 
For the (2, 1) system with L u = 1~ symmetry, these s v 
values correspond to K — 12, 12.3131 and 12.7, respec- 
tively. We find that there exists a value of the logarithmic 
derivative for each n that predicts energies as a function 
of Ro [via Eq. (|14l) ] that are in qualitative agreement with 
the dependence of the SV energies on r . For example, 
the SV energies of the lowest universal state for k = 12.7 
approach the ZR limit from above [pluses in Fig. Ob)], 
in agreement with the dependence of the energies ob- 
tained within the hyperspherical framework on Ro [see 
Fig. HoTd)]. For k = 12, in contrast, the SV energies of 
the ground state approach the ZR limit from below [cir- 
cles in Fig. HJa)], in agreement with the dependence of 
the energies obtained within the hyperspherical frame- 
work on Ro for small \L vq \ [see Fig.fTUTa)]. Although the 
comparison between the SV energies and those obtained 
within the hyperspherical framework leads to a consis- 
tent picture (including the existence or absence of bound 
states), the analysis unfortunately does not allow for the 
unambiguous extraction of the value of the logarithmic 
derivative or the generalized scattering length. 

To gain additional insight into the three-particle sys- 
tem, we analyze the structural properties of the (2, 1) 
system interacting through V g with L n = 1~, l/a s = 
and ro = 0.003a ho - Figure [TT] shows the scaled pair dis- 
tribution function P h i(r)r 2 for the heavy- light pairs for 
k = 1 through 12.5, while Fig. [T2l shows the scaled pair 
distribution function Phh{r)r 2 for the heavy- heavy pair 
for the same mass ratios. As k increases, the ampli- 
tudes of Ph\{T)r 2 and Phh(r)r 2 increase at small distances 
and decrease at large distances. The scaled pair distri- 
bution function for the spin-up — spin-down distance co- 
ordinate [75[ for k = 1 [dash-dotted lines in Fig. UTT a)]. 
e.g., shows a "two-bump structure" that has been previ- 
ously interpreted within an atom-dimer picture [60| : The 
"bump" at smaller r, r O.laho, reflects the formation 
of a dimer while the "bump" at larger r, r ~ 1.5aho, 
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FIG. 10: (Color online) Eigenenergy E of the trapped sys- 
tem at unitarity, obtained from Eq. (|14[) . as a function of Ro 
for (a) s v « 0.5579 (lowest allowed q value), (b) s„ = 1/2 
(lowest allowed q value), (c) s„ « 0.4180 (lowest allowed 
q value), and (d) s„ « 0.4180 (second lowest allowed g 
value). Dash-dash-dotted, dash-dotted, dotted, dashed and 
solid lines correspond to (a) L vq (xo) = 10, 1, 0, —1, and —10, 
(b) L vq {x ) = 2,1,0,-1, and -2, (c) L vq (x ) = 10,1,0,-1, 
and —2, and (d) L vq (xo) = 10,1,0,-1, and —10. 



reflects the fact that the "spare" spin-up atom sits fur- 
ther away from the spin-down atom than the spin-up 
atom that forms the dimer. For k = 6.7 [dotted lines 
in Figs. mT a) and lT27 a)]. which corresponds to the 40 K- 
6 Li mixture, the heavy-light and heavy-heavy pair distri- 
bution functions reflect the effective attraction between 
the two heavy atoms. Compared to the k = 1 system, 
the likelihood of finding the two like atoms at distances 
smaller than aho is significantly increased. For k — 12.5, 
both the scaled heavy-light and heavy-heavy pair distri- 
bution functions fall off approximately exponentially at 
length scales smaller than aho, as expected for a three- 
body bound state. 

Symbols in Fig. Q2] show the hyperradial density 




FIG. 11: (Color online) Dash-dotted, dotted, solid, and 
dashed lines show the scaled pair distribution function 
Phi(r)r 2 for the heavy- light pair of the (2, 1) system at uni- 
tary with L n = 1~ for (a) k — 1,6.7,11, and 11.5, and (b) 
K = 12, 12.3, 12.314, and 12.5, respectively. The heavy-light 



particles interact through V g (r) with ro = 0.003aho 
different scales of the axis in panels (a) and (b). 



Note the 



Phyper(R) for the FR interaction V g with r = 0.003ah o 
and X/a s = for various k. For comparison, lines show 
the hyperradial densities \F vq {R)\ 2 obtained from the ZR 
model. As can be seen from Eq. (T2"|) , F vq {R) depends on 
s v and q. For all k considered in Fig. 1131 we use the sq 
obtained by solving the hyperangular Schrodinger equa- 
tion for ZR interactions. The quantum number g is set 
to zero for k = 1 and 6.7 [see Eq. (fTTj) ]. For k > 8.619, 
we adjust the short-range boundary condition of F vq (R) 
so as to reproduce the FR energies, i.e., we calculate q 
according to Eq. (TIT)]) with E taken from the FR SV cal- 
culation. Figure [13] shows good agreement between the 
hyperradial densities obtained from the SV calculations 
and those obtained within the ZR model. For k = 12, 
e.g., the admixture of the irregular solution go q is clearly 
reflected in the hyperradial density: Using q = 0, i.e., 
using fo q only, results in a notably poorer description of 
the system (not shown). 

So far, we have considered systems with infinitely large 
two-body s-wave scattering length. To shed further light 
on the (2, 1) system in the vicinity of the three-body 
resonance, we also performed calculations for finite a s . 
For the free-space (2, 1) system with L n = 1~ and posi- 
tive a s , a single universal three-body bound state whose 
energy is a few times that of the two-body system has 
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FIG. 12: (Color online) Dash-dotted, dotted, solid, and 
dashed lines show the scaled pair distribution function 
Phh{r)r 2 for the heavy-heavy pair of the (2, 1) system at uni- 
tary with L n — 1~ for (a) k — 1,6.7,11, and 11.5, and (b) 
k = 12, 12.3, 12.314, and 12.5, respectively. The heavy-light 
particles interact through V g (r) with ro = 0.003ah o - Note the 
different scales of the axis in panels (a) and (b). 



been predicted to exist for k rs 8.173 to 12.917, and two 
universal bound states have been predicted to exist for 
k w 12.917 to 13.606 The three-body bound states 
discussed here for the trapped system are distinctly dif- 
ferent from these universal states, i.e., they do not ap- 
proach those discussed in Ref. [46| when uj approaches 
0. Symbols in Fig. [14] show the energies of the trapped 
three-body system as a function of l/a s for the inter- 
action potential V g with ro = 0.005ah o for three differ- 
ent mass ratios k, k = 12 (below the three-body reso- 
nance), k = 12.314 (at the three-body resonance), and 
k = 12.5 (above the three-body resonance). We note 
that the three-body energies depend strongly on ro. As 
expected, for negative (positive) a s , the three-body en- 
ergy lies above (below) the corresponding energy at uni- 
tarity. For comparison, a solid line shows the energy of 
the trapped two-body system with ZR interactions. For 
small \a s \ (a s > 0), the two-body energy varies approx- 
imately as 1/cig. To obtain a semi-quantitative descrip- 
tion of the three-body energies on the positive scattering 
length side, we write E/Eh Q = ci(a,h /a s — c^) C2 . Fitting 
our three-body energies for a,h /a s > 0.5, we find that 
the c-i coefficient changes from 1.65 over 1.55 to 1.52 as 
k changes from 12 over 12.314 to 12.5. The resulting fits 
are shown by dotted lines in Fig. 031 




0.5 1 1.5 

R/(2 1/2 a h0 ) 



FIG. 13: (Color online) Hyperradial density Phypcr(-R) of the 
(2, 1) system with L n = 1~ at unitary. Symbols show the re- 
sults from the SV calculation where the heavy-light particles 
interact through V s (r) with ro = 0.003aho for (a) k — 1 (cir- 
cles), k = 6.7 (squares), k = 11 (diamonds), and ft = 11.5 
(triangles), and (b) k = 12 (circles), k = 12.3 (squares), 
K — 12.314 (diamonds), and k = 12.5 (triangles). Lines show 
the hyperradial density |J*b w (as)| 2 , Eq. (|12|) . The so entering 
into Eq. (|12[) is taken from the ZR model, q is set to for 
K = 1 and 6.7, and calculated according to Eq. (|1U|) with 
E taken from the SV calculations for the other k. Note the 
different scales of the axis in panels (a) and (b). 



As discussed in Ref. [48|], we also considered systems 
consisting of two heavy and two light fermions interacting 
through Vg with l/a s = and various angular momenta. 
The addition of the light particle to the (2,1) system 
does not, to within our numerical resolution, lead to the 
appearance of a new resonance. Thus, adding a light 
particle leaves the system properties largely unchanged. 
The next subsection shows that adding a heavy fermion 
to the (2, 1) system does lead to the appearance of a new 
resonance. 



B. Four-body resonances 

Unlike the hyperangular Schrodinger equation for the 
(2, 1) system with ZR interactions, that for the (3, 1) 
system with ZR interactions is not analytically solu- 
ble. Thus, we employ the SV approach to solve the full 
Schrodinger equation for FR interactions, and then ana- 
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FIG. 14: (Color online) Circles, squares and diamonds show 
the SV energies for the (2, 1) system with L n = 1~ interacting 
through V g with ro = 0.005ai lo as a function of asho/fls f° r 
k — 12, 12.314 and 12.5, respectively. Dotted lines show fits 
to the SV energies (see text for details). For comparison, the 
solid line shows the energy for two trapped atoms interacting 
through V zr . 



TABLE I: Extrapolated ZR so values for small two- 
component Fermi gases at unitarity. The uncertainty of the 
so values is estimated to be in the last digit reported. The 
(2, 2) energies for k = 1 are taken from Ref. [72j and the (3, 2) 
energy for k = 1 is taken from Ref. [731 ]. 



K 


so(3,l) 


«o(4,l) 


so (2, 2) 


so(2,2) 


so(2,2) 


so (3, 2) 




1+ 


0" 


0+ 


1" 


2+ 


1" 


1 


4.08 


6.45 


2.509 


4.598 


3.418 


4.958 


2 


3.86 


6.15 


2.575 


4.357 




4.90 


4 


3.51 


5.68 


2.754 


3.997 


3.478 


4.85 


6 






2.886 


3.705 




4.76 


8 


2.79 




2.947 


3.430 


3.326 




10 






2.939 


3.138 


3.225 





lyze the eigenenergies and structural properties. In cer- 
tain cases, we "back out" so and q, i.e., we extract the 
lowest eigenvalue of the hyperangular Schrddinger equa- 
tion and the eigenvalues of the radial Schrddinger equa- 
tion. 

The energetically lowest lying state of the (3,1) sys- 
tem has L u = 1 + symmetry for the two-body scattering 
lengths a s of interest. In the a s — > 0~ limit, this can be 
verified by constructing the wave functions of the non- 
interacting system. We find numerically that this is also 
true at unitarity. Figures I15f a) and (b) show the en- 
ergy of the (3, 1) system as a function of ro for different 
k. Similarly to the (2, 1) system, the extrapolated ZR 
energy of the (3,1) system is approached approximately 
linearly from above for k = 1. For larger n (i.e., n = 4 
or 8), the ZR limit is approached approximately linearly 
from below. Table U and squares in Fig. \W\ summarize 
the extrapolated ZR energies for the (3,1) system with 
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FIG. 15: (Color online) Circles, pluses, squares and triangles 
show the SV energies for the (3, 1) system with L n = 1 + 
interacting through V s at unitarity as a function of ro for (a) 
K = 1,2,4 and 8, and (b) k = 10, 10.4, 10.5 and 10.6. Dotted 
lines are shown as a guide to the eye. 



L = 1 + symmetry for n < 8. The dependence of the 
energies on the range ro increases as K increases from 1 
or 2 to about 10, and then becomes comparatively small 
for k fa 10.5 [see squares in Fig. IToT b)]. For yet larger k, 
the energy varies roughly as 1/tq and approaches large 
negative values for small ro. In this regime, four-body 
bound states exist. Although our calculations for larger 
k, k > 8, do not allow for a reliable extrapolation of the 
energies to the ro — > limit, the existence of a four-body 
or (3, 1) resonance around k ps 10.4 or 10.5 is evident. 

Although our SV calculations for the (3,1) system are 
restricted to significantly larger ro than those for the 
(2, 1) system, it is clear that the behavior of the (3, 1) 
energies in the vicinity of the four-body resonance is qual- 
itatively similar to that of the (2, 1) energies in the vicin- 
ity of the three-body resonance. The main difference is 
that the four-body resonance occurs at a smaller mass ra- 
tio than the three-body resonance. Importantly, the SV 
energies for the Gaussian two-body potential V g for the 
(3, 1) system depend — as do those for the (2, 1) system — 
comparatively weakly on r in the vicinity of the reso- 
nance [squares in Fig. 115( b)]. The qualitatively similar 
dependence of the energies of the (3,1) and (2,1) sys- 
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FIG. 16: (Color online) Circles, squares and diamonds show 
the extapolated ZR energies as a function of k for the (2, 1) 
system with L n — V symmetry, the (3, 1) system with L n — 
1 + symmetry, and the (4, 1) system with L n = 0~ symmetry, 
respectively. As a guide to the eye, dotted lines connect data 
points. For the (3, 1) and (4, 1) systems, our FR SV energies 
allow for a reliable extrapolation to the ro — » limit only for 
relatively small n. 



terns on ro suggests that the four-body resonance occurs 
when so is approximately equal to 1/2. More specifi- 
cally, applying the framework detailed in Sec. Ill B[ i.e., 
assuming that the hyperradial and the hyperangular mo- 
tion separate not only for ZR interactions but also for 
FR interactions, the relatively weak dependence of the 
(3, 1) energies for k s» 10.5 implies that the four-body 
resonance occurs when sq is approximately 1/2. This 
interpretation assumes separability or approximate sepa- 
rability of the hyperangular and hyperradial parts of the 
wave function. If the coupling between different channel 
functions was appreciable for the FR interactions con- 
sidered in our SV calculations, then the hyperspherical 
framework of Sec. Ill Bl would need to be modified. While 
it is possible that the dependence of our SV energies on 
ro is consistent with a value of sq notably different from 
1/2 and non- vanishing channel coupling, we believe that 
the former scenario (i.e., approximate separability and 
so ~ 1/2 in the vicinity of the four-body resonance) is 
more likely. This conclusion is supported by our analysis 
of the hyperradial densities. 

Symbols in Fig. [17] show the hyperradial density 
-Phyper(-R) for the (3, 1) system with k = 1 — 10.6 in- 
teracting through the FR two-body potential V g with 
various ro. For k < 8, we find that the hyperradial den- 
sities for the Gaussian potential V g (symbols) are well 
reproduced by \F 0q (R)\ 2 , Eq. ifTZ]). with q set to zero 
and So determined by fitting the SV densities to Eq. (fl~2"|) 
[solid lines in Fig. flTT a)]. Alternatively, we set q = 
and use the SV energies to determine so via Eq. (fTT|l . 
For k = 1 and 4, the resulting hyperradial densities are 
indistinguishable from the lines shown Fig. ll~7T a). For 
k = 8, however, the description that treats sq as a fit- 
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FIG. 17: (Color online) Hyperradial density Ph ypC r(-R) of the 
(3, 1) system with L n = 1 + at unitary. Symbols show the 
results from the SV calculation where the heavy-light particles 
interact through V g (r) for (a) k = 1 and ro = 0.03ah o (circles), 
k = 4 and ro ~ 0.044ah o (diamonds), and k — 8 and ro = 
0.04aho (triangles); (b) k — 10 and ro ~ 0.054ai lo (circles), 
k — 10 and ro ~ 0.040ah o (diamonds), and k = 10 and ro ~ 
0.027ah o (triangles); and (c) k = 10.4 and ro ~ 0.041<2ho 
(circles), k — 10.5 and ro ~ 0.027ah o (diamonds), and k = 
10.6 and ro ~ 0.027ah o (triangles). In panel (a), lines show 
the hyperradial density \Fo q (R)\ 2 , Eq. (|12[) . with q set to zero 



and so determined by fitting the SV density to Eq. (|12p . 
the different scales of the axis in panels (a)-(c). 



Note 



ting parameter leads to a better description. Symbols 
in Fig. [TTTb) show the hyperradial densities for k = 10 
and three different ro values. As ro decreases, the the 
hyperradial density moves to larger R values and ap- 
proaches the hyperradial density expected for universal 
states. Symbols in Fig. fTTT c) show the hyperradial den- 
sities for k = 10.4, 10.5 and 10.6. As k increases, the 
maximum of Phypcr(-R) increases and moves to smaller R 
values. Moreover, the tail of Phypcr(-R) for K = 10.4—10.6 
starts to resemble that of a bound state. The ZR model, 
Eq. (Tl~2"|) , qualitatively but not quantitatively reproduces 
the FR hyperradial densities for n = 10 — 10.6 if the nor- 
malization constant No q , q and so are treated as fitting 
parameters. Based on a detailed analysis of the hyper- 
radial densities of the (2,1) system for various ro, we 
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believe that the ZR model would reproduce the FR hy- 
perradial densities of the (3, 1) system quantitatively if 
ro was smaller. 

Our calculations suggest that the sq for the (3,1) sys- 
tem with L n = 1 + symmetry changes notably over a 
small range of mass ratios. In particular, we find so ~ 1/2 
for k = 10.4 or 10.5. Recently, Castin and coworkers [76( 
reported that sq becomes zero for n = 13.384. While our 
so values are not in direct contradiction with this finding, 
an extrapolation of our results suggests that the Sq value 
of the (3, 1) system goes to zero at a smaller mass ratio 
than that found by Castin and coworkers. Future studies 
need to investigate this question in more detail. 

We also treated the energetically lowest-lying state of 
the (4, 1) system at unitarity, which has L n — 0~ symme- 
try. For this system, our calculations are restricted to the 
universal regime. Diamonds in Fig. [16] show the extrap- 
olated ZR energies and Table Q] shows the corresponding 
so values. We find that the hyperradial densities for the 
FR interaction potential (not shown) are well reproduced 
by Eq. (TT2")) with q = and So determined through a fit 
or by the SV energy. For completeness, Table U shows 
selected sq values for the (3, 1), (2, 2) and (3, 2) systems. 



IV. CONCLUSIONS 

We considered small fermionic two-component systems 
under external spherically symmetric confinement with 
1 < k < 13.607, i.e., in regimes where three-body Efimov 
physics is absent. Our calculations employed two differ- 
ent interaction potentials between the heavy and light 
atoms, a FR and a ZR interaction potential, and focused 
primarily on the unitary regime where the interspecies s- 
wave scattering length a s becomes infinitely large. Like 
atoms were assumed to be non-interacting. To address 
questions Q1-Q5 (see Sec. H]), we performed SV calcula- 
tions for few-fermion systems that interact through the 
FR interaction potential V s , and analyzed a subset of 
the eigenspectrum and selected structural properties such 
as the hyperradial density as a function of ro- The nu- 
merical results were interpreted within a hyperspherical 
framework that becomes exact in the ZR limit. We de- 
fined the generalized energy- dependent scattering length 
V S „(E), which characterizes the TV-body system in free 
space. This generalized scattering length is related to 
the hypervolume introduced in Ref. [77] to characterize 
bosonic three-particle systems and can be interpreted as 
being the result of an effective iV-body force that acts at 
small (or vanishing) hyperradii R. We used the gener- 
alized scattering length to connect the properties of the 
trapped V-body system with those of the free-space sys- 
tem. 

We found that the (2,1) and (3,1) systems at unitarity 
interacting through V g exhibit three-body and four-body 
resonances when k ~ 12.314 and k ~ 10.4, respectively. 
For the (2, 1) system with L n = l - symmetry, this mass 
ratio corresponds to, as has been shown by an indepen- 



dent calculation (4||, an so value of approximately 1/2. 
For the (3, 1) system with L n = 1 + symmetry, the So 
value is not known independently. We argued that the So 
value of the (3, 1) system is approximately 1/2 for around 
K = 10.4. It would be extremely valuable if this could be 
confirmed by an independent calculation that solves the 
hyperangular Schrodinger equation directly and does not 
rely, as our analysis, on backing so out from the full solu- 
tion. The three- and four-body resonances discussed here 
are obtained for the interaction potential V g . While the 
occurance of V-body resonances does depend, in general, 
on the details of the underlying two-body potential, we 
believe that the results presented here also apply to other 
short-range interaction potentials. In particular, we ar- 
gued that V-body resonances occur most likely when the 
so value that characterizes the V-body problem is ap- 
proximately 1/2. Experimentally, the value of so can 
possibly be tuned by varyin g th e lattice confinement or 
the intraspecies interactions [44{. While our calculations 
for the (4, 1) system were restricted to k values for which 
the system behaves universally, i.e., for which iV-body 
resonances are absent, recent work by Gandolfi and Carl- 
son [4^| found that the (4, 1) system in free-space exhibits 
an V-body resonance for n m 9.5 — 9.8. So far, V-body 
resonances have not been observed for larger systems and 
future research needs to address whether or not such res- 
onances exist. Our results for the (2, 1) and (3, 1) systems 
show, in agreement with Ref. [49| , that the energy of the 
iV-body bound state changes more rapidly for the (3, 1) 
system than for the (2,1) system for the same change Ak 
of the mass ratio. If iV-body resonances exist for larger 
systems, this trend is expected to continue. 

Our calculations show that FR effects become increas- 
ingly more important as the mass ratio k increases. The 
dependence on the range of the underlying two-body po- 
tential signals the breakdown of universality and has im- 
portant implications for theory and experiment. Since 
FR effects can be appreciable, theoretical treatments of 
unequal-mass Fermi systems based on ZR interactions 
may not be, in general, sufficient. Our calculations sug- 
gest that an analysis of experimental results for few- or 
many-fermion systems with sufficiently large mass ratio 
needs to account for non-universal FR effects. Near an 
iV-body resonance, the system's lifetime will be reduced 
due to the formation of molecules. If a two-component 
Fermi gas is prepared in a regime where ./V-body bound 
states exist, collapse is expected to set in. While this 
is, in many instances, an unwanted phenomenon, careful 
tuning in the vicinity of an .V-body resonance might open 
the possibility to study novel physics. On resonance, the 
./V-body clusters have a size comparable to that of s-wave 
dimers. At the (2, 1) resonance, e.g., the trimer is large. 
While bound tetramers and pentamers (and presumably 
larger clusters) exist, their wave functions have negligi- 
ble overlap with that of the trimer. This suggests that 
many-body systems with competing dimer and trimer in- 
teractions have finite lifetimes and may thus be prepared 
experimentally. Similarly, many-body systems with com- 
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peting two-body and four-body interactions should be ex- gratefully acknowledged, 
perimentally accessible. The study of these many-body 
systems is expected to uncover novel physics. 
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